Microbial isolation and characterization from two flex lines from the urine processor assembly onboard the international space station.

Urine, humidity condensate, and other sources of non-potable water are processed onboard the International Space Station (ISS) by the Water Recovery System (WRS) yielding potable water. While some means of microbial control are in place, including a phosphoric acid/hexavalent chromium urine pretreatment solution, many areas within the WRS are not available for routine microbial monitoring. Due to refurbishment needs, two flex lines from the Urine Processor Assembly (UPA) within the WRS were removed and returned to Earth. The water from within these lines, as well as flush water, was microbially evaluated. Culture and culture-independent analysis revealed the presence of Burkholderia, Paraburkholderia, and Leifsonia. Fungal culture also identified Fusarium and Lecythophora. Hybrid de novo genome analysis of the five distinct Burkholderia isolates identified them as B. contaminans, while the two Paraburkholderia isolates were identified as P. fungorum. Chromate-resistance gene clusters were identified through pangenomic analysis that differentiated these genomes from previously studied isolates recovered from the point-of-use potable water dispenser and/or current NCBI references, indicating that unique populations exist within distinct niches in the WRS. Beyond genomic analysis, fixed samples directly from the lines were imaged by environmental scanning electron microscopy, which detailed networks of fungal-bacterial biofilms. This is the first evidence of biofilm formation within flex lines from the UPA onboard the ISS. For all bacteria isolated, biofilm potential was further characterized, with the B. contaminans isolates demonstrating the most considerable biofilm formation. Moreover, the genomes of the B. contaminans revealed secondary metabolite gene clusters associated with quorum sensing, biofilm formation, antifungal compounds, and hemolysins. The potential production of these gene cluster metabolites was phenotypically evaluated through biofilm, bacterial-fungal interaction, and hemolytic assays. Collectively, these data identify the UPA flex lines as a unique ecological niche and novel area of biofilm growth within the WRS. Further investigation of these organisms and their resistance profiles will enable engineering controls directed toward biofilm prevention in future space station water systems.

Prediction of in-hospital mortality: An adaptive severity-of-illness score for a tertiary ICU in South Africa.

Background: A scoring system based on physiological conditions was developed in 1984 to assess the severity of illness. This version, and subsequent versions, were labelled Simplified Acute Physiology Scores (SAPS). Each extension addressed limitations in the earlier version, with the SAPS III model using a data-driven approach. However, the SAPS III model did not include data collected from the African continent, thereby limiting the generalisation of the results. Objectives: To propose a scoring system for assessing severity of illness at intensive care unit (ICU) admission and a model for prediction of in-hospital mortality, based on the severity of illness score. Methods: This is a prospective cohort study which included patients who were admitted to an ICU in a South African tertiary hospital in 2017. Logistic regression modelling was used to develop the proposed scoring system, and the proposed mortality prediction model. Results: The study included 829 patients. Less than a quarter of patients (21.35%; n=177) died during the study period. The proposed model exhibited good calibration and excellent discrimination. Conclusion: The proposed scoring system is able to assess severity of illness at ICU admission, while the proposed statistical model may be used in the prediction of in-hospital mortality. Contributions of the study: This study is the first to develop a model similar to the SAPS III model, based on data collected in South Africa. In addition, this study provides a potential starting point for the development of a model that can be used nationally.

Maternal iron status in early pregnancy and DNA methylation in offspring: an epigenome-wide meta-analysis.

BACKGROUND: Unbalanced iron homeostasis in pregnancy is associated with an increased risk of adverse birth and childhood health outcomes. DNA methylation has been suggested as a potential underlying mechanism linking environmental exposures such as micronutrient status during pregnancy with offspring health. We performed a meta-analysis on the association of maternal early-pregnancy serum ferritin concentrations, as a marker of body iron stores, and cord blood DNA methylation. We included 1286 mother-newborn pairs from two population-based prospective cohorts. Serum ferritin concentrations were measured in early pregnancy. DNA methylation was measured with the Infinium HumanMethylation450 BeadChip (Illumina). We examined epigenome-wide associations of maternal early-pregnancy serum ferritin and cord blood DNA methylation using robust linear regression analyses, with adjustment for confounders and performed fixed-effects meta-analyses. We additionally examined whether associations of any CpGs identified in cord blood persisted in the peripheral blood of older children and explored associations with other markers of maternal iron status. We also examined whether similar findings were present in the association of cord blood serum ferritin concentrations with cord blood DNA methylation. RESULTS: Maternal early-pregnancy serum ferritin concentrations were inversely associated with DNA methylation at two CpGs (cg02806645 and cg06322988) in PRR23A and one CpG (cg04468817) in PRSS22. Associations at two of these CpG sites persisted at each of the follow-up time points in childhood. Cord blood serum ferritin concentrations were not associated with cord blood DNA methylation levels at the three identified CpGs. CONCLUSION: Maternal early-pregnancy serum ferritin concentrations were associated with lower cord blood DNA methylation levels at three CpGs and these associations partly persisted in older children. Further studies are needed to uncover the role of these CpGs in the underlying mechanisms of the associations of maternal iron status and offspring health outcomes.

Long-stay medical-surgical intensive care unit patients in South Africa: Quality of life and mortality 1 year after discharge.

BACKGROUND: Although mortality is the primary measure of critical care outcome, the health-related quality of life (HRQOL) of survivors is often diminished. There is a paucity of South African research on HRQOL in intensive care unit (ICU) survivors. OBJECTIVES: To evaluate the 1-year post-discharge data of long-stay ICU patients, a group known to consume 20 - 40% of ICU resources. METHODS: A 1-year prospective observational study was conducted in a multidisciplinary medical-surgical ICU. Adult patients who were mechanically ventilated beyond 6 days were included. Clinical and mortality data were collected. Pre-admission and 6- and 12-month HRQOL were measured with the Short Form-36 questionnaire. Physical and mental component summary scores (PCS and MCS) were calculated. Associations between 12-month mortality and poor HRQOL scores were determined. RESULTS: Of 119 patients enrolled, 40.3% had sustained trauma, 19.3% were post-surgical and 40.3% had medical conditions; 29.2% were HIV-positive (HIV status was known for 74.8% of the cohort). The hospital and 12-month mortality rates were 42.9% and 57.4% (n=66/115), respectively. Age, longer ICU stay, higher disease severity scores and vasopressor use were associated with 12-month mortality. The survivors' median PCS and MCS at 6 and 12 months were significantly lower compared with pre-admission scores (both p<0.001). At 12 months, 53.1% of survivors demonstrated a poor PCS and 42.9% a poor MCS. Associations with poor 12-month PCS included longer ICU stay, male gender and trauma, while trauma and sepsis were associated with a poor 12-month MCS. Among the 19 trauma survivors, 78.9% had a poor MCS and/or PCS. Of previously employed patients, 54.8% were unemployed at 12 months. CONCLUSIONS: Patients ventilated beyond 6 days in a multidisciplinary ICU had a high mortality. Poor HRQOL at 12 months post discharge was frequently observed among survivors. Trauma was associated with poor 12-month outcomes. These findings highlight the need to further explore the outcomes of long-stay ICU patients in Africa.

Valproate-Associated Hyperammonaemic Encephalopathy.

We present the case of a 55 year old who presented multiple times with altered conscious levels. He was often treated as being post-ictal, when in fact, he had Sodium Valproate induced hyperammonaemic encephalopathy. Sodium Valproate can frequently increase ammonia levels, and in some patient lead to hyperammonaemic encephalopathy.

Ultra-sharp asymmetric Fano-like resonance spectrum on Si photonic platform.

In this paper, we report the generation of an ultra-sharp asymmetric resonance spectrum through Fano-like interference. This generation is accomplished by weakly coupling a high-quality factor (Q factor) Fabry-Pérot (FP) cavity and a low-Q factor FP cavity through evanescent waves. The high-Q FP cavity is formed by Sagnac loop mirrors, whilst the low-Q one is built by partially transmitting Sagnac loop reflectors. The working principle has been analytically established and numerically modelled by using temporal coupled-mode-theory (CMT), and verified using a prototype device fabricated on the 340 nm silicon-on-insulator (SOI) platform, patterned by deep ultraviolet (DUV) lithography. Pronounced asymmetric resonances with slopes up to 0.77 dB/pm have been successfully measured, which, to the best of our knowledge, is higher than the results reported in state-of-the-art devices in on-chip integrated Si photonic studies. The established theoretical analysis method can provide excellent design guidelines for devices with Fano-like resonances. The design principle can be applied to ultra-sensitive sensing, ultra-high extinction ratio switching, and more applications.

Maternal depression during pregnancy and cord blood DNA methylation: findings from the Avon Longitudinal Study of Parents and Children.

Up to 13% of women may experience symptoms of depression during pregnancy or in the postpartum period. Depression during pregnancy has been associated with an increased risk of adverse neurodevelopmental outcomes in the child and epigenetic mechanisms could be one of the biological pathways to explain this association. In 844 mother-child pairs from the Avon Longitudinal Study of Parents and Children, we carried out an epigenome-wide association study (EWAS) to investigate associations between prospectively collected data on maternal depression ascertained by the Edinburgh Postnatal Depression Scale in pregnancy and DNA methylation in the cord blood of newborn offspring. In individual site analysis, we identified two CpG sites associated with maternal depression in the middle part of pregnancy. In our regional analysis, we identified 39 differentially methylated regions (DMRs). Seven DMRs were associated with depression at any time point during pregnancy, 7 associated with depression in mid-pregnancy, 23 were associated with depression in late pregnancy, and 2 DMRs were associated with depression throughout pregnancy. Several of these map to genes associated with psychiatric disease and brain development. We attempted replication in The Generation R Study and could not replicate our results. Although our findings in ALSPAC suggest that maternal depression could be associated with cord blood DNA methylation the results should be viewed as preliminary and hypothesis generating until further replicated in a larger sample.

Proton therapy of prostate cancer by anterior-oblique beams: implications of setup and anatomy variations.

Proton therapy of prostate by anterior beams could offer an attractive option for treating patients with hip prosthesis and limiting the high-dose exposure to the rectum. We investigated the impact of setup and anatomy variations on the anterior-oblique (AO) proton plan dose, and strategies to manage these effects via range verification and adaptive delivery. Ten patients treated by bilateral (BL) passive-scattering proton therapy (79.2 Gy in 44 fractions) who underwent weekly verification CT scans were selected. Plans with AO beams were additionally created. To isolate the effect of daily variations, initial AO plans did not include range uncertainty margins. The use of fixed planning margins and adaptive range adjustments to manage these effects was investigated. For each case, the planned dose was recalculated on weekly CTs, and accumulated on the simulation CT using deformable registration to approximate the delivered dose. Planned and accumulated doses were compared for each scenario to quantify dose deviations induced by variations. The possibility of estimating the necessary range adjustments before each treatment was explored by simulating the procedure of a diode-based in vivo range verification technique, which would potentially be used clinically. The average planned rectum, penile bulb and femoral heads mean doses were smaller for initial AO compared to BL plans (by 8.3, 16.1 and 25.9 Gy, respectively). After considering interfractional variations in AO plans, the target coverage was substantially reduced. The maximum reduction of V 79.2/D 95/D mean/EUD for AO (without distal margins) (25.3%/10.7/1.6/4.9 Gy, respectively) was considerably larger than BL plans. The loss of coverage was mainly related to changes in water equivalent path length of the prostate after fiducial-based setup, caused by discrepancies in patient anterior surface and bony-anatomy alignment. Target coverage was recovered partially when using fixed planning margins, and fully when applying adaptive range adjustments. The accumulated organs-at-risk dose for AO beams after range adjustment demonstrated full sparing of femoral heads and superior sparing of penile bulb and rectum compared to the conventional BL cases. Our study indicates that using AO beams makes prostate treatment more susceptible to target underdose induced by interfractional variations. Adaptive range verification/adjustment may facilitate the use of anterior beam approaches, and ensure adequate target coverage in every fraction of the treatment.

Asymmetric split H-shape nanoantennas for molecular sensing.

In this paper we report on a very sensitive biosensor based on gold asymmetric nanoantennas that are capable of enhancing the molecular resonances of C-H bonds. The nanoantennas are arranged as arrays of asymmetric-split H-shape (ASH) structures, tuned to produce plasmonic resonances with reflectance double peaks within the mid-infrared vibrational resonances of C-H bonds for the assay of deposited films of the molecule 17ß-estradiol (E2), used as an analyte. Measurements and numerical simulations of the reflectance spectra have enabled an estimated enhancement factor on the order of 105 to be obtained for a thin film of E2 on the ASH array. A high sensitivity value of 2335 nm/RIU was achieved, together with a figure of merit of approximately 8. Our experimental results were corroborated using numerical simulations for the C-H stretch vibrational resonances from the analyte, superimposed on the plasmonic resonances of the ASH nanoantennas.

Investigation of time-resolved proton radiography using x-ray flat-panel imaging system.

Proton beam therapy benefits from the Bragg peak and delivers highly conformal dose distributions. However, the location of the end-of-range is subject to uncertainties related to the accuracy of the relative proton stopping power estimates and thereby the water-equivalent path length (WEPL) along the beam. To remedy the range uncertainty, an in vivo measurement of the WEPL through the patient, i.e. a proton-range radiograph, is highly desirable. Towards that goal, we have explored a novel method of proton radiography based on the time-resolved dose measured by a flat panel imager (FPI). A 226 MeV pencil beam and a custom-designed range modulator wheel (MW) were used to create a time-varying broad beam. The proton imaging technique used exploits this time dependency by looking at the dose rate at the imager as a function of time. This dose rate function (DRF) has a unique time-varying dose pattern at each depth of penetration. A relatively slow rotation of the MW (0.2 revolutions per second) and a fast image acquisition (30 frames per second, ~33 ms sampling) provided a sufficient temporal resolution for each DRF. Along with the high output of the CsI:Tl scintillator, imaging with pixel binning (2 × 2) generated high signal-to-noise data at a very low radiation dose (~0.1 cGy). Proton radiographs of a head phantom and a Gammex CT calibration phantom were taken with various configurations. The results of the phantom measurements show that the FPI can generate low noise and high spatial resolution proton radiographs. The WEPL values of the CT tissue surrogate inserts show that the measured relative stopping powers are accurate to ~2%. The panel did not show any noticeable radiation damage after the accumulative dose of approximately 3831 cGy. In summary, we have successfully demonstrated a highly practical method of generating proton radiography using an x-ray flat panel imager.

A WHO reference reagent to standardize haemagglutination testing for anti-A and anti-B in serum and plasma: international collaborative study to evaluate a candidate preparation.

BACKGROUND AND OBJECTIVES: The aim of the study was to evaluate a lyophilized serum preparation, 14/300, for its suitability to serve as a World Health Organization (WHO) Reference Reagent to standardize and control haemagglutination titrations for anti-A and anti-B in serum and plasma, in an international collaborative study. MATERIALS AND METHODS: Serum preparation 14/300 and two plasma-based reserve preparations, 14/304 (high titre anti-A) and 14/208 (high titre anti-B), were titrated by 24 laboratories in 13 countries using direct (DRT) and indirect (IAT) haemagglutination techniques. RESULTS: There was eightfold to 64-fold variation in reported titres per preparation and method across laboratories, that is, titres extended over 4-7 dilutions, although intralaboratory variability was generally good, with over 90% of replicate titres within a twofold range. There was a reduction in interlaboratory variability when titres of the reserve preparations were adjusted relative to those of the candidate Reference Reagent. CONCLUSION: The establishment of 14/300 as a WHO Reference Reagent for high titre anti-A and anti-B in serum, with nominal anti-A and anti-B titres of 128 for DRT, and nominal anti-A and anti-B titres of 256 for IAT, will facilitate global standardization of haemagglutination titrations for anti-A and anti-B in patient samples and blood components.

Stopping Criteria for Log-Domain Diffeomorphic Demons Registration: An Experimental Survey for Radiotherapy Application.

A crucial issue in deformable image registration is achieving a robust registration algorithm at a reasonable computational cost. Given the iterative nature of the optimization procedure an algorithm must automatically detect convergence, and stop the iterative process when most appropriate. This paper ranks the performances of three stopping criteria and six stopping value computation strategies for a Log-Domain Demons Deformable registration method simulating both a coarse and a fine registration. The analyzed stopping criteria are: (a) velocity field update magnitude, (b) mean squared error, and (c) harmonic energy. Each stoping condition is formulated so that the user defines a threshold ∊, which quantifies the residual error that is acceptable for the particular problem and calculation strategy. In this work, we did not aim at assigning a value to e, but to give insights in how to evaluate and to set the threshold on a given exit strategy in a very popular registration scheme. Experiments on phantom and patient data demonstrate that comparing the optimization metric minimum over the most recent three iterations with the minimum over the fourth to sixth most recent iterations can be an appropriate algorithm stopping strategy. The harmonic energy was found to provide best trade-off between robustness and speed of convergence for the analyzed registration method at coarse registration, but was outperformed by mean squared error when all the original pixel information is used. This suggests the need of developing mathematically sound new convergence criteria in which both image and vector field information could be used to detect the actual convergence, which could be especially useful when considering multi-resolution registrations. Further work should be also dedicated to study same strategies performances in other deformable registration methods and body districts.

Radiological Protection in Cone Beam Computed Tomography (CBCT). ICRP Publication 129.

The objective of this publication is to provide guidance on radiological protection in the new technology of cone beam computed tomography (CBCT). Publications 87 and 102 dealt with patient dose management in computed tomography (CT) and multi-detector CT. The new applications of CBCT and the associated radiological protection issues are substantially different from those of conventional CT. The perception that CBCT involves lower doses was only true in initial applications. CBCT is now used widely by specialists who have little or no training in radiological protection. This publication provides recommendations on radiation dose management directed at different stakeholders, and covers principles of radiological protection, training, and quality assurance aspects. Advice on appropriate use of CBCT needs to be made widely available. Advice on optimisation of protection when using CBCT equipment needs to be strengthened, particularly with respect to the use of newer features of the equipment. Manufacturers should standardise radiation dose displays on CBCT equipment to assist users in optimisation of protection and comparisons of performance. Additional challenges to radiological protection are introduced when CBCT-capable equipment is used for both fluoroscopy and tomography during the same procedure. Standardised methods need to be established for tracking and reporting of patient radiation doses from these procedures. The recommendations provided in this publication may evolve in the future as CBCT equipment and applications evolve. As with previous ICRP publications, the Commission hopes that imaging professionals, medical physicists, and manufacturers will use the guidelines and recommendations provided in this publication for implementation of the Commission's principle of optimisation of protection of patients and medical workers, with the objective of keeping exposures as low as reasonably achievable, taking into account economic and societal factors, and consistent with achieving the necessary medical outcomes.

Patient-specific stopping power calibration for proton therapy planning based on single-detector proton radiography.

A simple robust optimizer has been developed that can produce patient-specific calibration curves to convert x-ray computed tomography (CT) numbers to relative stopping powers (HU-RSPs) for proton therapy treatment planning. The difference between a digitally reconstructed radiograph water-equivalent path length (DRRWEPL) map through the x-ray CT dataset and a proton radiograph (set as the ground truth) is minimized by optimizing the HU-RSP calibration curve. The function of the optimizer is validated with synthetic datasets that contain no noise and its robustness is shown against CT noise. Application of the procedure is then demonstrated on a plastic and a real tissue phantom, with proton radiographs produced using a single detector. The mean errors using generic/optimized calibration curves between the DRRWEPL map and the proton radiograph were 1.8/0.4% for a plastic phantom and -2.1/ - 0.2% for a real tissue phantom. It was then demonstrated that these optimized calibration curves offer a better prediction of the water equivalent path length at a therapeutic depth. We believe that these promising results are suggestive that a single proton radiograph could be used to generate a patient-specific calibration curve as part of the current proton treatment planning workflow.

Validation of a deformable image registration technique for cone beam CT-based dose verification.

PURPOSE: As radiation therapy evolves toward more adaptive techniques, image guidance plays an increasingly important role, not only in patient setup but also in monitoring the delivered dose and adapting the treatment to patient changes. This study aimed to validate a method for evaluation of delivered intensity modulated radiotherapy (IMRT) dose based on multimodal deformable image registration (dir) for prostate treatments. METHODS: A pelvic phantom was scanned with CT and cone-beam computed tomography (CBCT). Both images were digitally deformed using two realistic patient-based deformation fields. The original CT was then registered to the deformed CBCT resulting in a secondary deformed CT. The registration quality was assessed as the ability of the dir method to recover the artificially induced deformations. The primary and secondary deformed CT images as well as vector fields were compared to evaluate the efficacy of the registration method and it's suitability to be used for dose calculation. plastimatch, a free and open source software was used for deformable image registration. A B-spline algorithm with optimized parameters was used to achieve the best registration quality. Geometric image evaluation was performed through voxel-based Hounsfield unit (HU) and vector field comparison. For dosimetric evaluation, IMRT treatment plans were created and optimized on the original CT image and recomputed on the two warped images to be compared. The dose volume histograms were compared for the warped structures that were identical in both warped images. This procedure was repeated for the phantom with full, half full, and empty bladder. RESULTS: The results indicated mean HU differences of up to 120 between registered and ground-truth deformed CT images. However, when the CBCT intensities were calibrated using a region of interest (ROI)-based calibration curve, these differences were reduced by up to 60%. Similarly, the mean differences in average vector field lengths decreased from 10.1 to 2.5 mm when CBCT was calibrated prior to registration. The results showed no dependence on the level of bladder filling. In comparison with the dose calculated on the primary deformed CT, differences in mean dose averaged over all organs were 0.2% and 3.9% for dose calculated on the secondary deformed CT with and without CBCT calibration, respectively, and 0.5% for dose calculated directly on the calibrated CBCT, for the full-bladder scenario. Gamma analysis for the distance to agreement of 2 mm and 2% of prescribed dose indicated a pass rate of 100% for both cases involving calibrated CBCT and on average 86% without CBCT calibration. CONCLUSIONS: Using deformable registration on the planning CT images to evaluate the IMRT dose based on daily CBCTs was found feasible. The proposed method will provide an accurate dose distribution using planning CT and pretreatment CBCT data, avoiding the additional uncertainties introduced by CBCT inhomogeneity and artifacts. This is a necessary initial step toward future image-guided adaptive radiotherapy of the prostate.

Validation of automatic contour propagation for 4D treatment planning using multiple metrics.

The aim of this work is to provide insights into multiple metrics clinical validation of deformable image registration and contour propagation methods in 4D lung radiotherapy planning. The following indices were analyzed and compared: Volume Difference (VD), Dice Similarity Coefficient (DSC), Positive Predictive Value (PPV) and Surface Distances (SD). The analysis was performed on three patient datasets, using as reference a ground-truth volume generated by means of Simultaneous Truth And Performance Level Estimation (STAPLE) algorithm from the outlines of five experts. Significant discrepancies in the quality assessment provided by the different metrics in all the examined cases were found. Metrics sensitivity was more evident in presence of image artifacts and particularly for tubular anatomical structures, such as esophagus or spinal cord. Volume Differences did not account for position and DSC exhibited criticalities due to its intrinsic symmetry (i.e. over- and under-estimation of the reference contours cannot be discriminated) and dependency on the total volume of the structure. PPV analysis showed more robust performance, as each voxel concurs to the classification of the propagation, but was not able to detect inclusion of propagated and ground-truth volumes. Mesh distances could interpret the actual shape of the structures, but might report higher mismatches in case of large local differences in the contour surfaces. According to our study, the combination of VD and SD for the validation of contour propagation algorithms in 4D could provide the necessary failure detection accuracy.

Computer models to study uterine activation at labour.

Improving our understanding of the initiation of labour is a major aim of modern obstetric research, in order to better diagnose and treat pregnant women in which the process occurs abnormally. In particular, increased knowledge will help us identify the mechanisms responsible for preterm labour, the single biggest cause of neonatal morbidity and mortality. Attempts to improve our understanding of the initiation of labour have been restricted by the inaccessibility of gestational tissues to study during pregnancy and at labour, and by the lack of fully informative animal models. However, computer modelling provides an exciting new approach to overcome these restrictions and offers new insights into uterine activation during term and preterm labour. Such models could be used to test hypotheses about drugs to treat or prevent preterm labour. With further development, an effective computer model could be used by healthcare practitioners to develop personalized medicine for patients on a pregnancy-by-pregnancy basis. Very promising work is already underway to build computer models of the physiology of uterine activation and contraction. These models aim to predict changes and patterns in uterine electrical excitation during term labour. There have been far fewer attempts to build computer models of the molecular pathways driving uterine activation and there is certainly scope for further work in this area. The integration of computer models of the physiological and molecular mechanisms that initiate labour will be particularly useful.

Interplay effects in proton scanning for lung: a 4D Monte Carlo study assessing the impact of tumor and beam delivery parameters.

Relative motion between a tumor and a scanning proton beam results in a degradation of the dose distribution (interplay effect). This study investigates the relationship between beam scanning parameters and the interplay effect, with the goal of finding parameters that minimize interplay. 4D Monte Carlo simulations of pencil beam scanning proton therapy treatments were performed using the 4DCT geometry of five lung cancer patients of varying tumor size (50.4-167.1 cc) and motion amplitude (2.9-30.1 mm). Treatments were planned assuming delivery in 35 × 2.5 Gy(RBE) fractions. The spot size, time to change the beam energy (τes), time required for magnet settling (τss), initial breathing phase, spot spacing, scanning direction, scanning speed, beam current and patient breathing period were varied for each of the five patients. Simulations were performed for a single fraction and an approximation of conventional fractionation. For the patients considered, the interplay effect could not be predicted using the superior-inferior motion amplitude alone. Larger spot sizes (σ ~ 9-16 mm) were less susceptible to interplay, giving an equivalent uniform dose (EUD) of 99.0 ± 4.4% (1 standard deviation) in a single fraction compared to 86.1 ± 13.1% for smaller spots (σ ~ 2-4 mm). The smaller spot sizes gave EUD values as low as 65.3% of the prescription dose in a single fraction. Reducing the spot spacing improved the target dose homogeneity. The initial breathing phase can have a significant effect on the interplay, particularly for shorter delivery times. No clear benefit was evident when scanning either parallel or perpendicular to the predominant axis of motion. Longer breathing periods decreased the EUD. In general, longer delivery times led to lower interplay effects. Conventional fractionation showed significant improvement in terms of interplay, giving a EUD of at least 84.7% and 100.0% of the prescription dose for the small and larger spot sizes respectively. The interplay effect is highly patient specific, depending on the motion amplitude, tumor location and the delivery parameters. Large degradations of the dose distribution in a single fraction were observed, but improved significantly using conventional fractionation.